Comparative study of neutron and nuclear matter 
with simplified Argonne nucleon-nucleon potentials 
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We present calculations of the energy per particle of pure neutron and symmetric nuclear matter with 
simplified Argonne nucleon-nucleon potentials for different many-body theories. We compare critically the 
Brueckner-Hartree-Fock results to other formalisms, such as the Brueckner-Bethe-Goldstone expansion up to 
third order, Self-Consistent Green's Functions, Auxiliary Field Diffusion Monte Carlo, and Fermi Hyper Netted 
Chain. We evaluate the importance of spin-orbit and tensor correlations in the equation of state and find these 
to be important in a wide range of densities. 
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I. INTRODUCTION 

The properties of homogeneous nuclear and neutron matter 
at high density play a crucial role in the determination of the 
structure of neutron star interiors QJJ] . Terrestrial nuclei pro- 
vide little input to constrain the equation of state (EOS) un- 
der the extreme conditions of density and isospin asymmetry 
within neutron stars. A potentially safe way to obtain the EOS 
thus relies on microscopic many-body calculations based on 
realistic nucleon-nucleon (NN) interactions. Over the years, 
several many-body approaches have been developed to de- 
scribe neutron star (and nuclear) interiors. The different ap- 
proaches might have very different physical foundations, but 
their final result is generally the same: a prediction for the 
density dependence of the energy per particle in neutron or 
symmetric nuclear matter. Here we want to progress in our 
understanding of the EOS by comparing quantitatively the re- 
sults provided by different many-body approaches. Similar 
benchmark calculations have taken place within the few-body 
community and have provided vital insight into the approxi- 
mation schemes at play yfl. 

The full operatorial structure of current high-quality NN 
potentials is presently too sophisticated for some state-of-the- 
art many-body schemes. Simplified versions of these poten- 
tials are therefore useful for benchmarking purposes. In par- 
ticular, we will use a family of simpler versions of the widely 
used Argonne Vis potential [3]. The V 8 ', VL and V[ potentials 
101 are built by removing operatorial components of the inter- 
action, while the remaining terms are readjusted (as indicated 
by the prime) in order to preserve as many lowest-order phase 
shifts and deuteron properties as possible. These potentials 
have been used in calculations of both infinite matter and nu- 
clei JH]. A natural question then arises: how well can these 
truncated potentials replace the original V^? 

To clarify this and other issues, we will first analyze in de- 
tail the properties of the family of Argonne potentials by in- 
specting their phase shifts in different partial waves. We will 
also examine the deuteron properties as predicted by these in- 



teractions. We will then study the EOS within the Brueckner- 
Hartree-Fock (BHF) many-body approach [6], which can han- 
dle straightforwardly all the variants of the Argonne potential. 
The BHF results are particularly insightful, because the to- 
tal energy can be directly connected to the partial wave ex- 
pansion and therefore to the microscopic properties of the 
in-medium NN interaction. We will also compare, when 
possible, the BHF results with other many-body approaches 
O, 0]> namely the Brueckner-Bethe-Goldstone approach up 
to third order in the hole-line expansion (BBG) ML the Self- 
Consistent Green's Function method (SCGF) IU0I - IT2I] . the 
Auxiliary Field Diffusion Monte Carlo (AFDMC) H], the 
Green's Function Monte Carlo (GFMC) [14], and the Fermi 
Hyper Netted Chain (FHNC) [15]. The comparisons should 
be helpful in quantifying theoretical uncertainties with respect 
to the EOS. As we shall see, the symmetric nuclear matter pre- 
dictions are particularly susceptible to the missing spin-orbit 
correlations in simplified potentials. 

This work complements and extends previous investiga- 
tions of some of the authors lfl6l Wfa . Our basic conclusions 
are related to the spin-orbit and tensor components of the NN 
interaction and NN in-medium correlations. Let us stress, in 
addition, that we do not attempt to obtain a "realistic" descrip- 
tion of the neutron or symmetric nuclear matter EOS. Our goal 
is solely to compare constructively the results obtained with 
different two-body potentials and many-body methods. A de- 
tailed calculation of the nuclear EOS would need to consider 
additional effects, particularly three-body forces 111811 . which 
are beyond the scope of the present work. 

The paper is organized as follows. In Sees. II and III we 
show the phase shifts and the deuteron properties predicted by 
the different Argonne potentials. Section IV is devoted to the 
BHF results of nuclear matter, both symmetric matter (SM) 
and neutron matter (NM), whereas in Sec. V we compare the 
results obtained with the different many-body approaches. Fi- 
nally, our conclusions are presented in Sec. VI. 
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II. PHASE SHIFTS 

The strong interaction part of the Argonne Vis potential can 
be expressed as a sum of 18 operators, 

Vij= £ Vpirijpfj. (1) 

p=l,18 

The first 14 operators are associated to the spin, isospin, ten- 
sor, spin-orbit, and quadratic spin-orbit components of the nu- 
clear force: 

P=W4 = l T .. T . 5 a a ( o, ■ G j ) I tj ■ X J ) . 
L-S,L-S(Ti-tj), 

L 2 , L 2 (Xi ■ Xj), L 2 (Oi ■ Oj), L 2 {Gj ■ Oj){Xi ■ Xj), 
(L-S) 2 ,(L-S) 2 (x r Xj). (2) 

The four additional operators, 

C 15 ^"" 18 = T 'h M°i ■ Oj), T^, (z Zl + r zj ) , (3) 

where 7y = 3t, ; t^ — T, • Xj is the iso-tensor operator, break 
charge independence. The radial functions that multiply each 
operator are adjusted by fitting experimental data on two-body 
scattering phase-shifts as well as deuteron properties. 

A family of simplified Argonne NN interactions has been 
devised to quantify the evolution of nuclear spectra with in- 
creasingly sophisticated NN interactions [4]. A given simpli- 
fied version, V„', is constructed by: i) eliminating the operato- 
rial structure with p > n and ii) refitting the remaining radial 
functions to reproduce as many properties of the original in- 
teraction as possible. Thus, V 8 includes all operators up to the 
spin-orbit term, but misses the components proportional to L 2 , 
(L ■ S) 2 , Tij, and % u + T zj . Similarly, V 6 ' is a NN force without 
spin-orbit (or further, p > 6) couplings. Finally, V4 does not 
even have a tensor coupling, but has been refitted to reproduce 
the binding energy of the deuteron. 

Realistic (or microscopic) NN interactions should fulfill a 
minimum set of requirements. In particular, realistic poten- 
tials are built to reproduce the Nijmegen database lfl9tl (which 
contains a full set of NN elastic scattering phase shifts up to 
energies of about 350 MeV) with an accuracy of % 2 /N^u ~ 1 ■ 
Only potentials that fulfill this condition should be used as in- 
put to the so-called ab initio many-body schemes, that aim at 
providing a first principles description of the EOS of neutron 
and symmetric matter. While Argonne Vis is, by all means, 
a realistic interaction, the Argonne family of simpler versions 
will necessarily violate this condition and thus they will be- 
come, in some sense, increasingly "less realistic". In partic- 
ular, one expects the reproduction of the phase shifts of high 
partial waves to be deteriorated as the operatorial structure of 
the interaction is simplified. One of our aims is to explore 
to which extent this deterioration has an impact on the EOS. 
This is particularly important in view of the fact that some 
approaches, at present, are limited to simplified forms of the 
Argonne family of potentials. In particular, the last generation 
of FHNC results has been computed with interactions up to V 8 



for both SM and NM. For ADFMC, the EOS of NM (SM) is 
only available with Vg (Vg). Diagrammatic approaches, such 
as BHF, BBG, or SCGF, have less limitations with respect to 
the structure of the original NN interactions. We will therefore 
provide results for Argonne Vis for these approaches. Within 
a given many-body approximation, calculations with V^, VL 
Vg, and V 4 ' should provide an indication of the importance of 
the missing operatorial components for the in-medium prop- 
erties. 

Fig. Q] shows the phase shifts of the lowest partial waves 
given by the different Argonne potentials J4J]. Vig, Vg, and 
Vg agree, by construction, in the 1 Sq, 3 Si, Pi, and 3 D[ par- 
tial waves. The over-simplistic Vj potential, however, does 
not have a tensor coupling, thus yielding zero mixing angles 
£i and £2. The 3 Di phase-shift is also badly reproduced for 
this potential, with an opposite sign relative to the other po- 
tentials. This is a direct consequence of the readjustment of 
the potential to bind the deuteron with the S channel only |0] . 
In the 3 £>i wave, the readjustment is such that it changes the 
nature of the interaction from repulsive (negative phase-shift) 
to attractive (positive phase-shift). 

For L — 1 , further discrepancies between phase shifts ap- 
pear. All potentials reproduce the l P\ phase-shift correctly, as 
the spin-orbit or tensor components are not active. Substantial 
differences, however, show up already in the 3 Po.i.2 waves. In 
these channels, the Vg and V 4 ' potentials deviate from the ex- 
perimentally fitted Vig results. In particular, the important 3 P2 
wave is grossly misrepresented with the Vg potential. In fact, 
this potential provides no correction in the S = 1 , T = 1 partial 
waves for the missing spin-orbit components and can hardly 
be considered realistic for this reason. 

Similarly large discrepancies are observed for the phase 
shifts of the l D2 and 3 £>2 partial waves. In the S = chan- 
nel, Vg and V 4 are identical and provide a too large phase-shift 
compared to Vis and Vg (which lie on top of each other in 
the plot). In the 5=1 channel, visible differences show up 
also between V 8 and Vig. For the L = 3 phase shifts, we have 
chosen to show the 3 f2 channel, where, again, substantial and 
visible differences appear between V 4 ', V 6 ', V 8 ', and Vis. 

The lowest right panels of Fig. Q] show the mixing parame- 
ters £1 and £2 of the i Si~ 3 D\ and 3 iV 3 f2 partial waves, respec- 
tively. On the one hand, the V 8 ' interaction reproduces very 
well the behavior of both parameters as given by the full Vi8, 
in spite of the slightly different 3 f2 phase shift. On the other 
hand, Vg' is only able to account for the low energy behavior 
(Slab _ 50 MeV) of £1 and £2. V 4 ' lacks a tensor component 
and therefore the mixing parameters associated to it are zero. 

Let us note, for further reference, that the BHF calculations 
are performed with partial waves up to J = 8. The SCGF re- 
sults presented here have been obtained for up to / = 4 in the 
T-matrix and up to J = 8 in the Hartree-Fock self-energy. For 
the BBG 3-hole calculations, all contributions have J < 5 and 
the convergence has been carefully tested. The differences 
in the J > 1 partial waves that we have just highlighted will 
therefore have an impact on the EOS. As a matter of fact, even 
phase-shift-equivalent potentials might predict different EOS 
due to their different (and physically unconstrained) off-shell 
structures. Naturally, one would expect these differences to 
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FIG. 1. (Color online) NN phase shifts for different potentials as a function of the energy in the laboratory E\ a ^. The solid lines represent the 
reference Argonne Vis results, whereas dotted, dashed, and dashed-dotted lines correspond to Vg, Vg, and V4, respectively. 



be small at low densities, where the physics is mainly domi- 
nated by L = components. As the density increases, how- 
ever, the differences in the higher partial waves start to show 
up. In non-perturbative diagrammatic calculations, the effects 
of high partial waves can be fed back to low momenta due to 
the self-consistency procedure. Likewise, relatively small dif- 
ferences in the phase shifts can have a substantial impact on 
the predictions of the EOS at densities even close to satura- 
tion. 

IIL DEUTERON PROPERTIES 

Even though the deuteron only explores the NN potential 
in the 3 5i- 3 Z)i partial waves, the analysis of the contributions 
of the different waves and operatorial components of the in- 
teraction provides a useful insight into the structure of the in- 
teraction Il20ll . We summarize the information related to the 
deuteron in Table [I] The first column gives the D-state proba- 
bility computed with different NN interactions. Although Po 
is not an observable, it provides an indication of the relative 
importance of the tensor coupling of the potential. By con- 



struction, the D-state probability of V g ' is the same as that of 
Vis, Pd = 5.78 %. For VL the probability decreases by less than 
10 %. As expected, Po = for V4, whose deuteron is a pure 
5-wave state. We also include the results of an additional po- 
tential, Vg, obtained by removing, without any readjustment, 
the spin-orbit components from Vg. For Vg, Pd is reduced by 
almost 20 %, indicating the importance of the spin-orbit com- 
ponents [p = 7,8 in Eq. ©J for the 3 Di wave and thus the 
ground state of the deuteron. Similar statements hold for the 
quadrupole moment of the deuteron, which is an observable. 
As such, Q c i is reproduced, by construction, with the original 
Vis force and the refitted V 8 ' and Vi forces. The lack of tensor 
components in V 4 ', however, implies a zero value of Q c /. The 
non-refitted interaction yields a value which is about 10 % 
larger than the experimental one. 

In columns 3, 4, and 5, we report the binding energy of 
the deuteron and its decomposition into kinetic and poten- 
tial terms. All the refitted potentials reproduce by construc- 
tion the total binding energy, E = —2.24 MeV. Note that this 
is the binding energy obtained only with the strong interac- 
tion components of the potential, i.e., when the small electro- 
magnetic terms are omitted. These repulsive electromagnetic 
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TABLE I. Deuteron £>-state probability Pp, quadrupole moment Qj (in fm 2 ), total binding energy, kinetic and potential energy, and their 
decomposition in partial waves, for different potentials. All energies are given in MeV. 



Force 


Pd (%) 


Qd 


E 


T 


V 


Ts 


T D 


Vs 


V D 


2V SD 


v« 


5.78 


0.27 


-2.24 


19.86 


-22.10 


11.30 


8.56 


-3.95 


0.77 


-18.91 


vl 


5.78 


0.27 


-2.24 


19.86 


-22.10 


11.30 


8.56 


-3.95 


0.77 


-18.91 


K 


5.33 


0.27 


-2.24 


18.70 


-20.94 


11.38 


7.32 


-4.68 


1.38 


-17.64 


yi 


0.00 


0.00 


-2.24 


11.41 


-13.65 


11.41 


0.0 


-13.65 


0.0 


0.0 


v 6 


4.64 


0.30 


-1.46 


14.96 


-16.42 


9.10 


5.86 


-3.43 


1.14 


-14.14 



terms shift the binding energy to the true experimental value 
of E = -2.22 MeV 0]. It is also relevant to note that the 
charge-dependent terms of Vis (p = 15, . . . , 18), described in 
terms of an iso-tensor operator, have no contribution in the 
iso-singlet deuteron state. 

It is well known that the deuteron binding energy results 
from a cancellation between a large positive kinetic and a large 
negative potential energy. For Vis these amount to T = 19.86 
MeV and V = —22.10 MeV, respectively. By construction, 
V 8 ' reproduces the same values as Vis. For V 6 ' there is a small 
variation in T and V due to the fact that the 3 S\ and 3 Z?i par- 
tial waves are not exactly identical to those of the V\% and 
Vg potentials. In particular, the kinetic energy decreases by 
about 1 MeV. Vj is also able to reproduce the total binding 
of the deuteron, but with much smaller kinetic (and therefore 
less negative potential) energies. In contrast, the potential V(, 
looses binding energy and also produces very noticeable dif- 
ferences for the kinetic and potential energies. In other words, 
a straightforward elimination of the spin-orbit components, 
without further readjustments, has large effects for the binding 
energy. 

It is also illustrative to separate the contributions of the 
3 5i and 3 ZJi states to the total kinetic and potential ener- 
gies. Assuming that the deuteron is a properly normalized 
combination of the i Si and 3 £>i partial waves, we define 
the contributions of the S and D states to the kinetic energy, 
T s = { 3 Si\T\ 3 Si) and T D = ( 3 £>i|r| 3 Z>i), and to the poten- 
tial energy, V s = ( 3 Si|V| 3 Si) and V D = ( 3 £>i|V| 3 £>i). The 
latter also receives a contribution from the 3 Si- 3 Di mixing, 
Vsd = ( 3 5'i I V J 3 X> 1 ) . These contributions are listed, for the dif- 
ferent potentials, in columns 6 to 10 of TableQ] 

For Vis, Vg, and V 6 ', the largest contribution to the potential 
energy actually comes from the mixing term, Vsd- This ac- 
counts for more than 85% of the final value of the potential 
energy. As already mentioned above, for Vj the deuteron is a 
pure 5-wave state, and therefore To, Vd, and Vsd vanish and 
the binding energy is obtained by accumulating a lot of attrac- 
tion in Vs- In spite of the fact that the spin-orbit components 
of Vis and Vg act explicitly only in the 3 Dj partial wave, when 
these are eliminated without readjustments in V5, all the con- 
tributions to the binding energy (and not only Vd) are altered. 
This is due to the fact that, due to the tensor coupling, the 
deuteron is obtained non-perturbatively from a combination 
of S, D, and mixing matrix elements. One can therefore say 
that when the spin-orbit component is not taken into account, 
a large change is induced in the wave function of the deuteron. 
In a wider picture, these results illustrate how the elimination 



TABLE II. Contribution of different components of the potential to 
the binding energy of the deuteron. All energies are given in MeV. 
See text for details. 



Force 


Central 


Tensor 


Spin-orbit 


L 1 


Vis 


-4.45 


-16.62 


-3.75 


2.72 


v{ 


-4.45 


-16.62 


-1.02 


0.00 


K 


-5.25 


-15.69 


0.00 


0.00 


v{ 


-13.65 


0.00 


0.00 


0.00 


v 6 


-3.84 


-12.58 


0.00 


0.00 



of operatorial components can have a relatively large impact 
in the binding energy of nuclear systems. The changes and 
differences induced by such elimination become more appar- 
ent when different channels are analyzed separately. While 
the refitting procedure in the V„' family of potentials seems to 
cure most deficiencies in the case of the deuteron, no such 
analogous procedure has been implemented in infinite matter. 
Consequently, one expects that even refitted potentials have a 
significant influence on the EOS of the infinite system. 

Alternatively, one can obtain a quantitative estimation of 
the different components of the potential by examining their 
expectation values in the ground-state wave function of the 
deuteron. We have grouped the 18 components in 4 different 
sets: the first four operators (p = 1, ... ,4), the tensor com- 
ponents Sjj (p = 5,6), the spin-orbit components L S and 
(L ■ S) 2 (p = 7, 8 and p — 13, 14), and the quadratic orbital an- 
gular momentum components L 2 (p = 9,. .., 12). The group 
of charge-dependent terms, p = 15, . . . , 18, does not contribute 
to the deuteron, as explained above. For Vis, the results of 
this decomposition are presented in the first row of Table HI] 
As expected, the largest contribution corresponds to the ten- 
sor component. All contributions are attractive, except that 
proportional to L 2 , which is slightly repulsive. The spin-orbit 
contribution is non-negligible and amounts to 17 % of the total 
potential energy. 

The second row of TableHIlshows the results for Vg. The to- 
tal potential energy, i.e., the sum of all the components, is the 
same as for Vi8- However, since the contribution of the first six 
operators yield the same result as Vis and the quadratic orbital 
angular momentum contribution is zero, the contribution of 
the spin-orbit terms is reduced to — 1 .02 MeV. In other words, 
the absence of the repulsive L 2 component is compensated by 
a smaller attractive contribution of the spin-orbit components, 
which decreases to 4.6 % of the total potential energy. Since 
these spin-orbit effects are more important for Vis than for Vg, 
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any technical difficulties associated to the many-body treat- 
ment of the spin-orbit components in Vg might artificially ne- 
glect contributions which are important in the case of the real 
deuteron, coming from the full interaction. 

For Vg only the first six operators contribute (row 3 in Ta- 
ble |llj. These components are about 1 MeV different from 
those in Vig and produce a slightly smaller total potential en- 
ergy compared to the full results. As we have discussed pre- 
viously, for Vg there is an important loss of binding energy. 
It is interesting to notice the difference between the spin-orbit 
contribution provided by V g ' ( — 1 .02 MeV) and the expectation 
value of the L ■ S component of Vg in the ground-state wave 
function provided by V(,, which amounts to —0.58 MeV. Ob- 
viously, this difference is due to the different wave functions 
generated by the two potentials. The wave function associated 
to V(, does not contain the spin-orbit correlations that are gen- 
erated by solving the Schrodinger equation with V 8 '. There- 
fore, the evaluation of the expectation value of the L ■ S com- 
ponent of Vg in the wave function provided by V(, gives a poor 
estimate of the spin-orbit contribution of V 8 '. This suggests 
that it would be necessary to incorporate the spin-orbit corre- 
lations in the wave function to recover the full contribution. 
The deuteron results illustrate effectively the non-negligible 
significance of spin-orbit contributions to the binding energy. 



IV. NUCLEAR EOS IN THE BHF APPROACH 

We now proceed to compare the results for the EOS of SM 
and NM obtained within the BHF approach for the family of 
simplified Argonne potentials described above. The BHF ap- 
roach represents the lowest order within the BBG expansion 
[H. In this formalism, the ground-state energy of infinite 
matter is computed from a diagrammatic expansion, which is 
regrouped according to the number of independent hole-lines. 
Within the BHF approach, the energy is given by the sum of 
only two-hole-line diagrams, including the effect of two-body 
correlations through the in-medium two-body scattering G- 
matrix. This takes into account the effect of the Pauli princi- 
ple on the scattered particles and the in-medium potential felt 
by each nucleon. 

As shown in Refs. U, the contribution to the energy 
from three-hole-line diagrams (which account for the effect 
of three-body correlations) is minimized when the so-called 
continuous prescription 12111 for the single-particle potential 
is adopted. This is a strong indication of the convergence of 
the hole-line expansion and we have adopted this prescription 
in our BHF calculations. We would also like to mention that 
the G-matrix has been calculated both in r- and £-space with 
two independent numerical codes and an excellent agreement 
has been found at all densities. We also remind the reader that 
no three-body forces are included in our calculations. 

In general, due to their different off-shell behavior, phase- 
shift equivalent potentials might lead to different saturation 
curves of nuclear matter. Within BHF, the saturation points 
of phase-shift equivalent interactions lie on the so-called "Co- 
ester band," which provides an empirical correlation between 
the saturation energy and density B22I1 . Quite obviously, po- 
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FIG. 2. (Color online) Energy per particle for nuclear (upper panel) 
and neutron (lower panel) matter as a function of density calculated 
in the BHF approach with the different Argonne potentials. 



tentials with different phase-shifts will also predict different 
saturation properties for infinite matter. We shall see in the 
following that a "Coester-like" correlation holds, however, for 
the different members of the simplified Argonne family of in- 
teractions. 

In Fig. [2] we compare the total energies of SM (upper panel) 
and NM (lower panel) within the BHF approach for the dif- 
ferent Argonne NN forces [6]. Several striking features arise 
from these comparisons. To begin with, one immediately no- 
tices the relatively large differences between the Vig (circles) 
and the Vg (squares) results in both cases. This is rather sur- 
prising, especially in view of the close agreements between 
phase shift and deuteron results presented earlier. In both SM 
and NM, Vg predicts a much more attractive EOS compared to 
Vig. In contrast, the Vl results (diamonds) are fairly close to 
the Vig calculations, in spite of their relatively different phase 
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TABLE III. Energies per particle (in MeV) for SM and NM at differ- 
ent densities and for different interactions calculated within the BHF 
approximation. 





Vis 


K 


K 


K 


P (fm- J ) 


SM 


NM 


SM 


NM 


SM 


NM 


SM 


NM 


0.03 


-6.2 


5.0 


-6.3 


4.9 


-6.1 


5.3 


-7.4 


5.6 


0.05 


-9.0 


6.7 


-9.3 


6.6 


-8.8 


7.3 


-11.2 


7.8 


0.10 


-12.6 


9.7 


-13.7 


9.2 


-12.2 


11.2 


-17.7 


12.2 


0.17 


-15.7 


14.0 


-18.4 


12.4 


-15.4 


16.7 


-26.4 


18.2 


0.20 


-16.4 


16.1 


-20.0 


13.8 


-16.2 


19.2 


-29.7 


20.9 


0.30 


-16.0 


25.3 


-23.7 


19.4 


-17.0 


28.9 


-39.2 


31.5 


0.40 


-11.9 


37.8 


-25.5 


26.2 


-15.2 


40.5 


-46.5 


44.0 


0.50 


-4.6 


53.6 


-26.0 


33.9 


-11.5 


53.9 


-51.8 


58.4 



shifts and deuteron predictions. The V 4 ' results are relatively 
reasonable for NM, where the 3 5Di channel is absent, while 
they are clearly unrealistic for SM. As a matter of fact, sym- 
metric matter does not even saturate before p = 0.5 fm~ 3 for 
this potential. For a more detailed insight, we refer the reader 
to Table [HI] where the energies of SM and NM at different 
densities are listed for the Argonne family of potentials. 

Another relevant question associated to the spin-orbit cou- 
pling arises when comparing the results predicted by the Vg 
and % potentials. VL has been constructed specifically by 
modifying only the central potential in the T = 0, S = 1 states 
0]. As a consequence, both potentials are identical in the 
T = 1 partial waves and therefore they predict the same EOS 
for pure NM (see lower panel of Fig. |2j. The effect of cor- 
recting the central T = 0, S = 1 term is only seen in SM, and 
corresponds to a rather density-independent gain of binding of 
a few MeV with respect to the Vg results (see upper panel of 
Fig. 13. 

The BHF approach is almost unique within all the many- 
body approaches, in the sense that the total energy of the 
system can be linked to the partial wave expansion of the 
in-medium NN interaction. This provides an interesting in- 
sight into the impact of different partial waves on the EOS 
and has motivated us to use BHF as a reference calculation. 
An analysis of the different partial wave contributions to the 
total energy at p = 0.3 frrV 3 is listed in Table IIVI for SM and 
NM. In SM, the largest difference between the Vig and V 8 ' re- 
sults arises from the 3 SD\ wave (4.8 MeV), even though the 
phase shifts in these coupled channels are practically identical. 
We have verified that this effect is actually due to the differ- 
ent single-particle potentials. The differences in higher-order 
partial waves thus influence, via the single-particle states, the 
lowest partial waves. Other substantial contributions to the 
difference between Vis and V 8 ' stem from the 'So (1.3 MeV), 
the l D 2 (0.8 MeV), and the 3 D 2 (0.9 MeV) states. 

Somewhat surprisingly, the total energy of Vl is actually 
closer to that of Vis, even though the individual partial wave 
contributions are rather different. For instance, the 3 PF 2 en- 
ergy difference is very large (—12.3 MeV), in agreement with 
the unrealistic phase shifts of Vi in this channel. This large 
discrepancy, however, is cancelled by opposite (and also rel- 
atively large) differences in the 3 Pq and 3 Pi channel. We 



TABLE IV. Partial wave decomposition of the binding energy per 
particle (in MeV) of SM and NM at p = 0.3 fm~ 3 for different po- 
tentials. The total sum comprises partial waves up to J = 8. The row 
in brackets contains the partial sum of the 3 Pq, 3 Pi, 3 PF 2 states. 



SM 



State 


Vis 


vi 


K 


vl 


Vis 


vi 


vi 


vl 


T= 1 


















'So 


-21.9 


-23.2 


-22.4 


-22.1 


-20.4 


-22.4 


-20.8 


-20.5 


3 i>o 


-5.0 


-5.1 


-9.1 


0.6 


-4.2 


-4.4 


-11.4 


1.4 


3 P, 


20.2 


19.9 


15.6 


1.8 


31.5 


30.8 


22.6 


4.2 




-16.5 


-16.6 


-4.2 


3.1 


-25.4 


-25.7 


-3.6 


7.1 


(*P> 


-1.3 


-1.8 


2.3 


5.5 


1.9 


0.7 


7.6 


12.7) 


l D 2 


-5.9 


-6.7 


-6.7 


-6.7 


-10.0 


-12.1 


-12.1 


-12.0 


T = 


















3 SDi 


-20.4 


-25.2 


-23.9 


-43.5 










'Pi 


7.3 


7.2 


7.3 


7.3 










y D 2 


-8.0 


-8.9 


-8.8 


-4.8 










All 


-16.0 


-23.7 


-17.0 


-39.2 


25.3 


19.4 


28.9 


31.5 



NM 



specify the sum of all 3 P states in the fifth row of Table [TTT] 
This shows that, in the case of Vg, all 3 P waves amount to a 
difference of only —3.3 MeV in the total energy. Together 
with the opposite-sign differences in the 3 5Di (3.5 MeV) and 
3 £>2 channels (0.8 MeV), the total energy eventually becomes 
close to the Vis value. 

As we have just discussed, the binding energies (and satura- 
tion curves) of SM obtained with the Vl and the Vi potentials 
are rather different. These large differences are associated to 
the phase-shift non-equivalence of the two interactions. Note, 
in particular, that the largest differences between the two sets 
of results manifest themselves in the channels where the two 
potentials are actually not equivalent (i.e., where Vg has been 
modified with respect to V g '). The different mixing ofthe 3 S£>i 
channels induces a difference of 1.3 MeV. Once again, even 
though the 3 P individual contributions are rather different, the 
overall sum leads to a relatively small difference of 4.1 MeV. 
This is the largest contribution to the total binding energy dif- 
ference of 6.9 MeV. Since the energies associated to Vg are 
more attractive, a higher saturation density is found in com- 
parison with Vis and Vi. We conclude that an intricate series 
of compensations leads to closer agreement between Vg and 
Vis than between Vg and Vis. This is clearly misleading in 
view of the insufficient phase shift reproduction of Vg com- 
pared to Vg. 

A similar issue is observed in the partial wave decomposi- 
tion of the V4 results. The T = 1 partial waves disagree sub- 
stantially on a one-by-one basis (in some cases by more than 
10 MeV), but the overall sum of the 3 P waves differs from 
the Vis results by a smaller number, —6.8 MeV. The most ex- 
treme difference is obtained, as expected, in the 3 5Di channel, 



which is about 20 MeV more attractive for V4 than for the rest 
of potentials. Ultimately, it is this extreme additional bind- 
ing, due to the lack of partial-wave coupling and the extreme 
readjustment of the 3 Dj channel, that drives the non-saturating 
behavior of V4. 

The partial wave differences in NM (see left columns in 
Table |IVt also carry interesting information. On a channel- 
by-channel basis, the comparison between V^ and Vg results 
is more favourable than that between Vis and V 6 '. Yet, once 
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again, the discrepancies are reduced in the overall sum, so 
that the Vl final results are closer to Vig. At p = 0.3 fm -3 , 
NM is about 10 MeV more bound with the V 8 than with V 6 '. 
This difference is entirely due to the spin-orbit component. In 
particular, the phase-shift non-equivalence in the ^PF2 chan- 
nel is evidenced by a difference of more than 20 MeV. Sim- 
ilarly, the V4 results for NM are only slightly more repulsive 
than the others. The partial wave decomposition suggests that 
this is due to a reshuffling of partial wave contributions which 
would, separately, deviate substantially from more realistic re- 
sults. 

Finally, we would like to comment on the specific effect of 
the spin-orbit components in a neutron-rich medium. First, 
notice that in NM the differences between Vg' = Ve and V g ' are 
only due to the suppression of the spin-orbit components in 
Vg, as no other readjustments are applied. Clearly the elimi- 
nation (rather than the readjustment) procedure for spin-orbit 
components in the interaction does not have a small effect in 
the EOS of NM, as evidenced by the large differences between 
the Vg and V 8 ' results in Fig.|2] 

One might be tempted to attribute the small NM differences 
between the VL and Vig interaction to: i) small spin-orbit com- 
ponents in the NN interaction or ii) small spin-orbit correla- 
tions in the medium. As we have seen, however, the partial 
wave expansion suggests that this is a rather fortuitous coin- 
cidence. First, both potentials are not truly phase-shift equiv- 
alent. Second, the in-medium corrections are quite different, 
as evidenced by the different partial wave components of the 
energy. All in all, our findings suggest that spin-orbit com- 
ponents should not be arbitrarily eliminated in ab initio cal- 
culations. Extending the argument to V 8 , one could say the 
same for the remaining missing operators in the interaction. 
Their effect is relevant for the EOS and needs to be dealt with 
properly. 

A similar reasoning can be extended to SM. The results pro- 
vided by V(, (without readjusting the potential) and V 6 ' (with 
readjustment) are rather close in that case, indicating that the 
readjustment in the central T = 0, S = 1 channel is relatively 
small. One can therefore say that the differences between the 
SM V 6 ' and V 8 results are mainly due to the suppression of the 
spin-orbit component in the interaction itself and they give 
rise to rather large corrections. In other words, the suppres- 
sion of the spin-orbit component does not have a small effect 
on the EOS. The fact that the V[ 8 and Vg SM results are rel- 
atively close at low densities is only due to a cancelation ef- 
fect in different partial waves. This can hardly be ascribed to 
a smallness of spin-orbit interactions or in-medium correla- 
tions. Let us point also out that the situation for the EOS is 
different to that of light nuclei, where Vg, Vg, and Vis predict, 
in the framework of GFMC, somewhat similar spectra and to- 
tal energies [4]. However, it is rather difficult to find a simple 
correspondence between the EOS differences presented here 
and the GFMC for finite nuclei. 



V. COMPARISON WITH OTHER MANY-BODY METHODS 

Several approaches have been devised over the decades to 
treat the nuclear matter many -body problem. Critical compar- 
isons between the approaches can help us learn how correla- 
tions, produced by different pieces of the original NN poten- 
tials, are treated within each scheme 11231 12411 . The ultimate 
goal of these comparisons should not be to find a "good" or a 
"bad" EOS, but rather to understand and quantify the differ- 
ences among them. Eventually, a well-founded comparison 
might also delineate the theoretical uncertainties of present 
generation EOSs. 

The only way to carry out a meaningful comparison be- 
tween many-body approaches is by starting from the same 
underlying NN interaction. While in principle one would like 
these to be phase-shift equivalent and as realistic as possible, 
some many-body approaches are currently limited by the op- 
eratorial structure of the NN force. FHNC and AFDMC, in 
particular, are coined to treat Argonne-type interactions, in- 
volving the sum of products of local radial functions and op- 
erators. Some parts of the full Vi 8 interaction cannot yet be 
fully included in these many-body schemes. To be able to 
compare between many-body methods, we will perform cal- 
culations with different approaches with the same family of 
underlying NN potentials. By progressively adding compo- 
nents to the NN force, we also hope to elucidate the role of 
such operators in the different many-body treatments of the 
EOS. 

Results are summarized in Figures [3] and [4] which repre- 
sent the energy per particle as a function of density for both 
SM and NM, respectively. The four panels in each figure give 
results for different NN interactions. Within each panel, BHF 
results are shown as circles; SCGF as squares; FHNC as dia- 
monds; AFDMC as up triangles; and BBG as down triangles. 
In the case of NM we also show available GFMC (right tri- 
angles) results II 1411 for the Vg and V 8 potentials. Let us first 
comment the SM results and move later to NM results. 



A. Symmetric nuclear matter 

1. BBG 

We start by discussing the SM results obtained with the 
BBG approach at the three-hole line level. Three-hole line 
contributions are especially small when using the continuous 
prescription for the single-particle spectrum ]9|]. Using the 
gap prescription, the BBG value converges also to the same 
results |9l|24j]. In this case, however, the three-hole line con- 
tribution is sizeable and therefore the lowest-order BHF re- 
sults cannot be taken as a good estimation of the energy. We 
note here that the BBG results shown in this work have been 
calculated with the so-called /^-matrix ||6t] whereas, as it has 
been said before, the BHF ones have been obtained with the 
G-matrix. We have checked, however, that at the BHF level 
the differences between G-matrix and /f-matrix calculations 
are negligible for both SM and NM, and all the potentials. 
Only a difference of about 2 MeV is found in the case of the 
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FIG. 3. (Color online) Energy per particle of symmetric nuclear matter as a function of density calculated with several many-body approaches 
for different Argonne potentials: V4 (upper left), Vg (upper right), V 8 ' (lower left), and V\$ (lower right). 



Vg potential for both SM and NM at the highest density con- 
sidered. 

We confirm the good agreement between our continuous 
choice BHF and the BBG calculations for all the NN interac- 
tions considered. The results differ by a few Me V at large den- 
sity, indicating, as already mentioned, a rather good conver- 
gence of the hole-line expansion. We should mention, how- 
ever, that the differences between the BHF and BBG results 
are even smaller when the BHF is done with the /^-matrix 
instead of the G-matrix. Note that the three-hole line con- 
tribution can be either attractive or repulsive, depending on 
the interaction under consideration. In particular, it is min- 
imal for the most realistic potential, Vjg (lower right panel). 
The BBG predictions for Vg' represent minor repulsive cor- 
rections to the BHF results and both lie well below the other 
approaches. Note, in particular, that nuclear matter saturates 
beyond 0.5 fm~ 3 in the BHF case. As explained earlier, this is 
due to the partial wave contributions associated to spin-orbit 
and tensor forces, which are particularly strong for this poten- 
tial. This is confirmed by comparison with the V(, results in 
the same panel, in which these forces are removed and a very 
strong reduction of binding is observed. 

In spite of the self-consistency process required to calculate 
the energy, the BHF approach is not, strictly speaking, ther- 
modynamically consistent. The Hugenholtz-Van Hove theo- 
rem is typ ically violated by about 20 MeV in BHF calcula- 
tions ll25l l26ll . The latter theorem specifies that the chem- 
ical potential should equal the quasi-particle energy at Icf ■ 
Within the hole-line expansion, however, the single-particle 
spectrum is merely an auxiliary quantity and should therefore 



not be taken as a definition of the chemical potential. A way 
to deal with thermodynamical consistency within BHF theory 
has been discussed in Ref. 12711 . This should in principle be 
applicable to BBG calculations as well. 

2. SCGF 

Another way to approach the many-body problem is 
through the SCGF method [8]. In this case, a diagrammatic 
expansion is employed to solve for the in-medium one-body 
propagator, rather than for the energy of the system. For in- 
finite matter, the method is conventionally applied at the lad- 
der approximation level. With respect to the G-matrix, the 
in-medium ladder interaction presents two major differences. 
First, hole-hole intermediate states are considered in addi- 
tion to the typical BHF particle-particle propagators. Sec- 
ond, these intermediate propagators are fully dressed, i.e., 
expressed in terms of spectral functions rather than through 
single-particle energies only. 

At a formal level, the comparison between the BHF and the 
SCGF approaches is not straightforward. Even though both 
approaches arise from a diagrammatic expansion, the infinite 
subsets of diagrams considered in both approaches are not the 
same. Moreover, the summation procedures are also some- 
what different, with the Dyson equation being used in SCGF 
to dress all internal propagators. Whereas the BHF formalism 
in the continuous choice can be derived from the ladder SCGF 
formalism after a series of approximations 12811 . this is not the 
case for the full BBG expansion. 
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FIG. 4. (Color online) Same as Fig. [5] but for neutron matter. 



In principle, if both BBG and SCGF were carried out to 
all orders, they should yield identical results. BBG theory, 
however, is an expansion in powers of density (or hole-lines), 
and the three-hole line results seem to indicate that it con- 
verges quickly. The error in the SCGF expansion is more dif- 
ficult to quantify, as one cannot directly compute (or even es- 
timate) the values for the diagrams of other structural types. 
We show a few representative diagrams for the perturbation 
expansion of the energy in Fig. [3] Dashed (continuous) lines 
represent interactions (fermions). Up to second order, i.e., 
diagrams (a) and (b), both approaches include the same dia- 
grammatic contributions. The only discrepancies at this level 
would arise from the (potentially self-consistent) treatment of 
internal lines. Within the SCGF approach, this would auto- 
matically give rise to diagrams like (c), which, in BBG the- 
ory, are considered separately. At higher orders, ring diagram 
iterations, as in diagram (e) or (h), would be included in the 
three-hole line BBG expansion, but not in the SCGF ladder 
resummation. Similarly, diagram (i) represents a third-order 
bubble diagram that is not explicitly incorporated in the SCGF 
theory. Diagram (f) is a typical hole-hole scattering process 
included in SCGF and absent in BHF. 

It is well established that the differences between SCGF 
and BHF result in an overall repulsive effect to the binding 
energy of both SM |28l-[30ll and NM which is mainly 
due to the inclusion of hole-hole propagation. Phase space 
arguments suggest that this repulsive effect should increase 
with density. The dressing of intermediate propagators is cap- 
ital for the thermodynamical consistency of the method. The 
Hugenholtz-Van Hove theorem is well fulfilled in these calcu- 



lations. At the technical level, we would like to point out that 
the SCGF results presented here have been obtained as zero- 
temperature extrapolations of finite-temperature calculations. 
The details of the extrapolation procedure will be presented 
elsewhere ll3lll . Pure zero-temperature calculations have been 
obtained using different numerical techniques by the Ghent 
H29I1 and Krakow ll ill groups. 

Compared to BHF, the repulsive effects within the SCGF 
approach are quite sizable in SM, especially when the ten- 
sor component is considered. For the case of Vig (lower 
right panel), for instance, the saturation point shifts from 
[0.25 firT 3 ,- 16.8 MeV] for BHF to [0.17 firT 3 ,- 11.9 MeV] 
for SCGF. While the shift seems to go towards the right satu- 
ration density, the value of the SCGF saturation energy is quite 
high. A similar tendency (lower saturation density, quite less 
binding at saturation) is also found for Vg and V g ', in agree- 
ment with the findings of Ref. M2911 . The oversimplified V4 
does not seem to saturate within the SCGF approach either. 
Note that, for this potential, the results are very close to BHF, 
which could suggest that, if tensor and spin-orbit correlations 
are not present, the many-body problem is under better control 
(see also the upper left panel in Fig. |4j). 

In the SCGF approach, the total energy is usually obtained 
through the Koltun sum rule @]. Even though this gives direct 
access to the total true kinetic and potential energies, the sum 
rule cannot be explicitly analyzed in terms of partial waves, as 
we have done for BHF. It is therefore difficult to pin down the 
observed differences to specific terms in the NN interaction. 
Alternative ways of computing the energy, as those suggested 
in Ref. Ill 111 , could provide such a partial wave decomposition. 
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FIG. 5. Diagrams in the perturbation expansion of the total energy. Only direct terms are shown for simplicity. 



3. FHNC 

Several approaches based on the variational principle have 
also been devised to treat the nuclear matter problem i32tl . 
These are rather different from the methods based on non- 
perturbative sums of diagrams presented earlier, but the basic 
principle is rather simple. The expectation value of the Hamil- 
tonian, 



E = 



in a trial wave function, 



Pgr | H | gV) 




0F. 



ermi-Gas ; 



(4) 



(5) 



sym 



provides an upper bound to the exact energy. The trial wave 
function is built by means of a correlation operator, 



FiJ = I,f P (rtj)Otj 



(6) 



that describes the correlations induced by the NN interaction. 
In principle, this correlation operator should have the same 
operatorial structure as the interaction, see Eq. (fTJ. In practice, 
though, this is generally not the case, since the calculation of 
the expectation value for the full Hamiltonian is technically 
very involved. 

The total energy can be evaluated in a diagrammatic 
cluster expansion with the aid of the Fermi Hyper Netted 
Chain/Single Operator Chain (FHNC/SOC) integral equations 
ll23ll . which sum Meyer-type diagrams containing up to an in- 
finite number of nucleons. The sum is, however, incomplete, 
as some topologies and operatorial structures are difficult to 
include within infinite summations. For instance, the elemen- 
tary diagrams are generally not included in FNHC calcula- 
tions. Similarly, the spin-orbit correlations cannot be chained, 



and are usually evaluated at the three-body cluster level. On 
top of this variational upper bound to the energy, one can add 
second-order perturbative corrections (A£ CBF ) calculated in 
the framework of the correlated basis function (CBF) theory 



The FHNC results for SM with Vjj and V 6 ' are presented in 
the corresponding panels of Fig. [3] For our comparisons, we 
use the FHNC calculations reported in Refs. lfl5ll33ll . There 
is no clear tendency with respect to our BHF reference calcu- 
lations. Whereas the V 6 ' FHNC results are more attractive than 
the BHF calculations, a significantly repulsive behavior is ob- 
served for Vj>. Consequently, the saturation points are quite 
different to those of BHF. For V 6 ', the FHNC results saturate 
around [0.39 ftrT 3 ,-20.1 MeV], while the V g ' calculation sat- 
urates at [0.35 ftrr 3 , -17.6 MeV]. 

If the variational procedure had been exactly performed, the 
FHNC results should provide an upper bound to the energy 
per particle. However, as pointed out in Refs. lfl5ll . the un- 
certainties and difficulties in the treatment of the spin-orbit 
correlations in the FHNC formalism have hindered the inclu- 
sion of spin-orbit correlations in the wave function, i.e., the 
underlying NN interaction and correlation function are not 
treated on an equal footing. A widely debated question is if 
the variational character of FHNC results is preserved after the 
approximations involved in the calculational procedure. Re- 
cent comparisons between FHNC and AFDMC results have 
enlightened this basic question lll5ll . Notice also that the in- 
clusion of 2 particle-2 hole corrections calculated at second 
order in the framework of the CBF formalism would give an 
additional attraction to the FHNC results lfl7ll. 



Contrary to what we have observed within BHF, the FHNC 
results seem to be quite independent of the underlying inter- 
action. This has been attributed to the smallness of the spin- 
orbit components of the NN interaction as well as the cancel- 
lation of spin-orbit correlations in the in-medium wave func- 
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tion. These arguments stem from comparisons with AFDMC 
results [ 15]. The impossibility to switch off spin-orbit corre- 
lations in BHF hinders a direct comparison with the V 8 ' FHNC 
calculations, where the spin-orbit correlations have been omit- 
ted. As in the case of the deuteron, to have a quantitative 
idea of the contribution of the spin-orbit components, we re- 
port BHF results for the Vg interaction. The EOS of SM with 
Vg is shown in the lower left panel of Fig. [3] As explained 
earlier, the difference between the V 8 and Vg results is only 
caused by the spin-orbit components in the original interac- 
tion. The differences turn out to be rather large and therefore 
this comparsion, although not conclusive, does not support the 
idea that switching off spin-orbit correlations produces small 
changes in the energy. In any case, FHNC calculations with a 
more elaborate treatment of the spin-orbit correlations would 
be highly desirable. 



4. QMC 



series of limitations, particularly those associated to the sign 
problem. 

Let us finally comment on the results obtained with V 4 '. 
As an unrealistic interaction without a tensor component, V| 
yields an unrealistic EOS for SM which does not saturate for 
any of the many-body approaches discussed here. This po- 
tential, however, can be used for benchmarking purposes. As 
a matter of fact, in spite of its limitations, the V 4 ' results are 
useful, as we find a large degree of agreement among the 
different many-body approaches. This confirms the long ac- 
cepted notion that, in the nuclear matter case, the correlations 
induced by spin-orbit and tensor components of the interac- 
tion are probably on the basis of the large differences when 
comparing different many-body methods. As a matter of fact, 
the antisymmetrization procedure in NM eliminates part of 
these components (particularly, the tensor coupling) and, as 
we shall see in the following, this results into a much closer 
agreement among all the many-body calculations. 



Quantum Monte Carlo (QMC) methods are very success- 
ful in describing the ground state of infinite system of bosons, 
like atomic liquid 4 He HH, or fermions, like liquid 3 He i35ll . 
Additional efforts have allowed the extension of the QMC 
method to nuclear systems, which have more complicated 
interactions and correlation structures. However, the accu- 
racy of QMC in its different versions, be it AFDMC OJ] or 
GFMC [14], is limited by the fermion sign problem 13611 . Up 
to now, the safest way to deal with this problem in nuclear 
systems is to keep the sample walk within a fixed nodal sur- 
face. This is an approximation and, consequently, a potential 
limitation of the QMC approach. AFDMC and GFMC dif- 
fer in the way they treat the spin and isospin degrees of free- 
dom. AFDMC samples the spin-isospin states by introducing 
Hubbard-Stratonovich auxiliary fields, whereas GFMC sums 
them completely. This fact limits in a drastic way the number 
of nucleons that GFMC can consider, generally up to about 
16. The auxiliary field strategy allows AFDMC to efficiently 
sample spin-isospin correlations in systems with a sufficient 
number of nucleons (N = 114) to ensure that the finite size 
corrections are small. A recent comparison has demonstrated 
that both methods give very close results for neutron drops 
with AT < 16|37ll. 

In spite of its recent progress, the AFDMC method is still 
not able to work with the full Vjg potential. Technical prob- 
lems, again mainly associated with the spin-orbit structure of 
the interaction and the trial wave function, are the reason for 
the lack of AFDMC results for SM with potentials containing 
spin-orbit components (V^ or V 8 ). We will use the most recent 
version of the AFDMC results in our comparisons [38]. V 4 re- 
sults are rather similar to our BHF reference calculation up to 
0.2 fm~ 3 , but the AFDMC predictions become more attractive 
as density increases. The AFDMC results for SM with VL as 
shown in the upper right panel of Fig. [3] are quite close to 
the FHNC predictions. Note for instance that AFDMC with 
Vg saturates around [0.37 far 3 , -18. 8 MeV]. However, the 
agreement between both approaches should not be taken as a 
final consistency check. Note that these results still present a 



B. Neutron matter 

The NM results for different NN interactions and many- 
body approaches are presented in Fig. [4] Overall, it is fair 
to say that the differences between the different potentials are 
significantly smaller for NM than for SM. As just mentioned, 
we attribute this similarity to the fact that, for a fully isospin- 
polarized system, the number of active partial waves is re- 
duced. This is particularly true for the 3 SDi coupled wave, 
which is inactive in this system. Consequently, the differences 
between potentials are further reduced. Moreover, the lack of 
tensor correlations presumably facilitates the solution of the 
many-body problem. 

For V 4 ' (upper left panel), all reported approaches give quite 
similar results. This supports the idea that, when the potential 
is just central, a good agreement between the different many- 
body techniques is found. Let us in particular stress the fact 
that, up to p = 0.3 fm~ 3 , the results of all the different ap- 
proaches fall within a rather narrow window of less than 5 
MeV. 

When the tensor component of the force is taken into ac- 
count, as in Vg (see upper right panel), the differences up to 
p = 0.3fm~ 3 are still relatively small. At this density, we 
find e BHF = 28.6 MeV, e BBG = 26.9 MeV, e SCGF = 31.2 MeV, 
e afdmc = 25.5 MeV, and e FHNC = 25.0 MeV. In other words, 
the inclusion of the tensor component increases the uncer- 
tainty (measured as the spread predicted by these calculations) 
in the many-body calculations by 1 — 2 MeV. Note that, how- 
ever, a broader discrepancy is observed at high densities, with 
FHNC lying below all the other approaches. In general, SCGF 
provides the more repulsive results. As mentioned earlier, Vg 
and Vg produce the same results for neutron matter because 
the readjustment of the potential does not affect the NM par- 
tial waves. 

Since the results for Vg and Vg are the same for NM, the 
differences between V 8 and V 6 ' give a direct measure of the 
spin-orbit contributions in the NM channels where it is active. 
The incorporation of the spin-orbit components in V 8 ' produces 
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an overall attraction with respect to the V 6 ' results, apart from 
FHNC. The most attractive variation corresponds to BHF, 
while the SCGF, FHNC, and AFDMC approaches stay rather 
close to their Vg counterparts at all densities. Note that also 
GFMC results seem to have a stronger dependence on the 
underlying interaction, and that they are in good agreement 
with the BBG ones. At p = 0.3 fm~ 3 , we find e BHF = 18.3 
MeV, e BBG = 20.4 MeV, <? SCGF = 24.8 MeV, e AFDMC = 24.7 
MeV, and g FHNC = 26.4 MeV. The spread in these results 
has increased substantially compared to the simpler poten- 
tials. Again, the differences increase with density. Notice, 
however, that in this case the FHNC results stay above all the 
other methods for all densities, whereas for Vg this approach 
provided the most attractive results. 

We would like to stress once again that the spin-orbit contri- 
bution in the NN interaction produces a sizeable contribution 
for the NM BHF EOS with respect to Vg. This is in contrast to 
the small effect observed in the AFDMC or the FHNC meth- 
ods. For these two methods one should notice that spin-orbit 
correlations (and also the tensor ones in the AFDMC case) are 
not included in the trial wave function, although the AFDMC 
allows for their fully generation in the diffusion process. 

The results for Vis are only available for the three per- 
turbative diagrammatic approaches, BHF, SCGF, and BBG. 
The agreement between them is in general rather good. At 
0.3 fm- 3 , for instance, we find e BHF = 25.4 MeV, e BBG = 25.2 
MeV, and e SCGF = 29.4 MeV This represents a relatively 
small spread in values compared with the differences we have 
observed in SM. One should keep in mind that these differ- 
ences increase with density. Therefore, as a general remark we 
can conclude that in NM the change of the EOS with respect 
to the BHF results for both SCGF and BBG is less important 
than in SM. This is particularly true for the V' A and Vg po- 
tentials. For the V g ' potential, larger differences are observed 
between these three approaches even in NM. In all cases, the 
SCGF method produces the more repulsive EOS. 

VI. CONCLUSIONS 

By performing calculations of the energy per particle for 
nuclear and neutron matter within different many-body ap- 
proaches, we have investigated the properties of the equation 
of state as obtained with a family of frequently used Argonne 
potentials. The many-body approaches we have used are: the 
Brueckner-Hartree-Fock, the Brueckner-Bethe-Goldstone ap- 
proach up to third order in the hole-line expansion, the Self- 
Consistent Green's Function method, the Auxiliary Field Dif- 
fusion Monte Carlo, and the Fermi Hyper Netted Chain. We 
have analyzed critically the origin of underlying uncertainties 
in these calculations. The subtraction and refitting procedure 
of different operatorial structures, inherent in the family of Ar- 
gonne NN interactions used in the present paper, has provided 
us with substantial insight. 

At the phase shift and deuteron levels, we find that already 
the Vg potential produces fairly large deviations relative to the 



original V\% model. While at the two-body level Vg is almost 
a clone of Vis for S and P partial waves, we have found that, 
when included in many-body calculations, it produces rela- 
tively large differences. These discrepancies are driven by the 
phase-shift differences in higher partial waves as well as by 
off-shell effects. The Vg NN interaction should therefore be 
regarded critically in high-precision applications. 

The overall infinite-matter binding energies obtained with 
the Vg' interaction are actually closer to the Vis results, in spite 
of the unsatisfactory reproduction of the P-wave phase shifts. 
We believe that this agreement is, however, due to an artificial 
reshuffling of different partial wave contributions. To support 
this claim, we have presented the partial wave decomposition 
of the total energy in BHF calculations, which has indeed con- 
firmed an overall cancellation of larger differences in the total 
energy. In the T = 1 channel, this potential is not readjusted 
for the missing spin-orbit component, which otherwise pro- 
duces fairly important contributions to the binding energy. It 
might well be that microscopic properties, other than the EOS, 
are also affected by these large discrepancies. 

Finally, while the V 4 ' model is clearly unrealistic and should 
only be used for benchmarking purposes, we have found that 
it is at this level that the different many-body approaches agree 
more. This confirms the long accepted hypothesis that the ten- 
sor (and spin-orbit) components of the interaction and their 
in-medium treatment are at the heart of most of the observed 
discrepancies. Overall, one needs to consider these with more 
attention before drawing more definite conclusions. Studies 
like the one presented here help us in quantifying the uncer- 
tainties in state-of-the-art predictions of the EOS, originating 
either from the microscopic interaction or the many -body the- 
ory. Only a thorough understanding of these uncertainties will 
allow us to provide a well-founded connection between the 
physics of neutron stars, the fundamental strong interaction 
acting on its constituents, and the many-body correlations at 
play. 
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